MOD-^ CONVERGENCE 

FREDDY DELBAEN, EMMANUEL KOWALSKI, AND ASHKAN NIKEGHBALI 

Abstract. Using Fourier analysis, we study local limit theorems in weak-convergence prob- 
lems. Among many applications, we discuss random matrix theory, some probabilistic mod- 
els in number theory, the winding number of complex brownian motion and the classical 
situation of the central limit theorem, and a conjecture concerning the distribution of values 
f^ ■ of the Riemann zeta function on the critical line. 



'^ ■ 1. Introduction 



^ 



In p3] , the notion of mod-gaussian convergence was introduced: intuitively, it corresponds 

P^ I to a sequence of random variables X„ that - through the Fourier lens - "look like" a sum 

Oh ■ Gn + Yn where (G^) is a sequence of gaussian variables with arbitrary variance and (F„) 

^ ■ is a convergent sequence independent from G„. However, most interest lies in cases where 

c^ \ this simple-minded decomposition does not exist: what remains is the existence of a limiting 

function $, not necessarily a Fourier transform of a probability measure, such that the limit 

theorem 

^ : (1) lim E[e**^"]-iE[e**^"] = ^{t) 

^2 ! holds, locally uniformly, for t G M. 

\f^ • In the same spirit, we introduce in this paper a notion of "convergence" where the ref- 

[~>^ . erence law is not necessarily gaussian but a fairly general probability law, with integrable 

O \ characteristic function yj. Under suitable conditions, we are able to prove a general local 

limit theorem which extends the result recently found in [THl Th. 4]. 

As illustrations of the consequences of this framework, we mention two results which, to 

the best of our knowledge, are new: 



Theorem 1 (Local limit theorem for the winding number of complex brownian motion). 
For M > 0, let 9u denote the argument or winding number of a complex brownian motion W^ 
such that Wq = 1. Then for any real numbers a < b, we have 

lim -^ F[a<e^<b] = -{b-a). 

It— S>oo 2 TT 

This is proved in Section 13.21 
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Theorem 2 (Local limit theorem for unitary matrices). For n > 1, let Qn denote a random 
matrix which is Haar- distributed in the unitary group U{n). Then for any bounded Borel 
subset B G C with boundary of Lebesgue measure 0, and for any b & C, we have 



lim ^^^ P 



lo£C n 
\ogdetil-gn)-^beB 



2tt 



where m{-) denotes the Lebesgue measure on C. 



This, together with similar facts for the other families of classical compact groups, is 
proved in Section 13.71 

We emphasize that these two applications are just examples; this paper contains quite 
a few more, and it seems certain that many more interesting convergence theorems can be 
proved or understood using the methods of this paper. 

Notation and preliminaries. Our random variables will take values in R'^, a fixed 
(i— dimensional space, and we denote by |t| the Euclidian norm in R''. We will use the 
Landau and Vinogradov notations f = 0{g) and f <^ g in some places; these are equivalent 
statements, and mean that there exists a constant c > such that 



■'5 



|/(a;)| <cg{x) 

for all a; in a set X which is indicated. Any suitable value of c is called "an implied constant" , 
and it may depend on further parameters. 

A sequence of probability measures (/in)„ on R"' converges weakly to a probability mea- 
sure if for all bounded continuous functions / : M'^ — > R, we have lim„ / / (i/i„ = j f dfi. 
Equivalently we can ask that the convergence holds for C°° functions with compact support. 
Levy's theorem asserts that this is equivalent to the pointwise convergence of the charac- 
teristic functions Jexp{it ■ x) dfin — ^ fexp{it ■ x) dfi. Levy's theorem can be phrased as 
follows. If (fn is the sequence of characteristic functions of probability measures //„, if fn(t) 
converges pointwise to a function ip, if this convergence is continuous at the point 0, then 
(y9 is a characteristic function of a probability measure fi, fj,n converges weakly to /i and the 
convergence of (/?„ to (p is uniform on compact sets of R*^. We recall that the convergence is 
continuous if x„ — )■ in M'^ implies V5„(a;„) — )• ip{0) = 1. 

We say that a sequence of random variables converges in law if the image measures (or 
laws) converge weakly. Most of the time one needs a scaling of the sequence. This is 
for instance the case in the central limit theorem, which in an elementary form says that 
for a sequence of independent identically distributed real-valued random variables, (X„)n, 
E[X„] = 0,E[X^] = 1, the normalised (or rescaled) sequence ^^^'^^ converges weakly to 
the standard gaussian law. 

In the applications below we will use different kinds of scaling. In the higher- dimensional 
case, we will scale the random variables using a sequence of linear isomorphisms (or non- 
degenerate matrices) A^ : R'^ — )■ R'^. The inverse of these matrices will be denoted by 
Hn'- R'^ — )■ R'^. The transpose of a linear map or matrix A is denoted by A*. 

Our methods are based on Fourier analysis and we will use basic facts from this theory 
freely. We define the Fourier transform as is usually done in probability theory, namely 

f{t)= / exp{it ■ x)f{x)m{dx). 
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The inversion formula is, at least when / G //^(R*^), given by 



^xrf 



fi^) = hr / exp{-tt-x)fit)dt. 

\2,TX J J^d 

In particular, when /i is a probability measure with an integrable characteristic function 
(y9, we get that yU is absolutely continuous with respect to Lebesgue measure m, ft <t^ m, and 
its density is given by 

^(a;) = f ^ 1 / exp(-zt ■ x)(f{t) m{dt), 

which is therefore continuous. 

The proof of our main result (TheoremH]) is based on the following approximation theorem: 

Theorem 3. Suppose f : M.'^ ^ M. is a continuous function with compact support. Then for 
each rj > we can find two functions gi,g2'- M'^ — )■ M such that 

(1) ^,^ have compact support, 

(2) g2<f<gi, 

(3) f^d{gi-g2)it)m{dt) <ri. 

This is a standard result in Fourier analysis, see e.g. Bretagnolle and Dacunha-Castelle [5]; 
for the sake of completeness, we sketch a proof in Appendix A. 

2. MOD-Lp CONVERGENCE 

2.1. Definition. We now explain our generalization of the definition in [13]. First of all, 
we fix t/ > 1 and a probability measure fi onM.'^. We then assume given a sequence (X„) of 
random variables defined on a probability space (fi, J-", P) and taking values also in M*^. We 
define (fn to be the characteristic function of X„. We now consider the following properties: 

• HI. The characteristic function ip of the probability measure /i is integrable; in 
particular, /i has a density dfi/dm, with respect to Lebesgue measure m. 

• H2. There exists a sequence of linear automorphisms An G GL(i(M), with inverses 
S„ = A~^, such that S„ converges to and (/9„(S*t) converges continuously at 
(or what is equivalent: uniformly on compact sets) to v^(t). In other words, the 
renormalized random variables S„(X„) converge in law to /x. 

• H3. For all k > 0, the sequence 

fn,k = V^n(S*t)l|2*i|<A; 

is uniformly integrable on M'^; since fn,k are uniformly bounded in L^ and L°° (for 
fixed k), this is equivalent to the statement that, for all fc > 0, we have 

(2) lim sup/ \ipn{T,lt)\l\j:*^t\<km{dt) = 0. 

Remark 1. Property HI excludes discrete probability laws, such as Poisson random variables. 
However, similar ideas do apply for such cases. We refer to [17] (for the case of Poisson 
distributions) and to [2] (for much more general discrete distributions, where earlier work of 
Hwang [lOj are also relevant) for these developments. 
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Remark 2. Property H3 will typically be established by proving an estimate of the type 

(3) \ipn{Kt)\<h{t) 

for all n > 1 and all t eW^ such that |SJjt| < A;, where /i > is an integrable function on W^ 
(which may depend on k). 

We give a name to sequences with these properties: 

Definition 1 (Mod-*/? convergence). If /x is a probability measure on W^ with characteristic 
function yj, X„ is a sequence of M'^-valued random variables with characteristic functions yj^, 
and if the properties HI, H2, H3 hold, then we say that there is mod—if convergence for 
the sequence X„. 

Below, we will comment further on the hypotheses, and in particular give equivalent 
formulations of H3. In Section I3.H we also explain the relation with conditions arising in 
classical convergence theorems. 

To make the link with the original definition in |13], i.e., the assumption that a limit 
formula like ([1]) holds, we observe that mod-y? convergence will hold when HI is true and 
we have 

• H2'. There exists a sequence of linear automorphisms A^ G GLci(]R), with inverses 
S„ = A~^, such that E„ converges to 0, and there exists a continuous function 
$ : M'^ -)- C such that 

(4) ^„(t) = <|.(tMA;t)(l + o(l)) 

uniformly for t such that |S*t| < A;, for arbitrary A; > 0. 

In many applications considered in this paper (not all), this stronger condition holds, or 
is expected to hold. It is very likely that, when this is the case, the "limiting function" $ 
also carries significant information, as discussed already in special cases in [121 §4]. 

2.2. Local limit theorem. We now state and prove our main result, which is a local limit 
theorem that shows that, when mod— y? convergence holds, the expectations E[/(X„)] (for 
reasonable functions /) do not converge, but are well-controlled: they behave like 

|det(AOr'^(0) / fdm 
dm. Jf,d 

as n goes to infinity. The proof turns out to be very simple: 

Theorem 4 (Local limit theorem for mod-ip convergence). Suppose that mod—ip convergence 
holds for the sequence X„. Then we have 



det(A„)|E[/(X„)]^^(0) / fdm, 

dm J^d 



for all continuous functions with compact support. Consequently we also have 

(5) |det(A„)|P[X„Gfi]^^(0)m(5). 

for relatively compact Borel sets i? C M^ with m{dB) = 0, or in other words for bounded 
Jordan-measurable sets 5 C M'^. 



Proof. We first assume that / is such that / has compact support; using Theorem [3l the 
case of a general continuous function with compact support will follow easily. We write 



E[f{Xn)] = [ f{x)dfin{x) 



where yU„ is the law of X„, before applying the Plancherel formula and the inversion formula 
to transform this into 

By the linear change of variable t = S* s, we get 

E[/(X„)] = (27r)-Vet(S„)| / cp^Ks) HKs) m{ds) . 

Now fix k so that the support of / is contained in the ball of radius k; we then have 
E[/(XJ] = (27r)-Vet(S„)| [ MKs) KKs) m{ds) . 

J|S*s|<fc 

The integrand converges pointwise to (p{s)f{0) according to the assumption H2. The 
condition H3 of uniform integrability then implies the convergence in L^. One can see this 
quickly in this case: for any e > 0, and for any a > large enough, we have 

/ \^niKs)'i-\j:*„s\<kf iKs)\ m{ds) < \\f Woo / \(pniKs)'i-\i:*„s\<k\m{ds) < e 

^|s|>a J\s\>a 

for all n by ([2]). On \s\ < a, the pointwise convergence is dominated by ||/||ool|s|<a, hence 
/ </'n(S;s)l|s*.|<fc/(S;s)m(rfs) ^/(O) / ip{s)m{ds). 

J\s\<a J\s\<a 

For a large enough, this is /(O) J (f, up to error e, hence we get the convergence 



cpniKs)fiKs)m{ds)-^f{0) / ip{s)m{ds). 

|S*s|<fc JR-i 



-dj^^r.^ I ,./.N„./J„N ^/^ 



Finally, this leads to 

|det(A„)|E[/(X„)]^(27r)-V(0) / ^{s)m{ds) = ^{0) I f{s)m{ds), 

jRd am j^d 

which concludes the proof for / with / with compact support. 

Now if / is continuous with compact support, we use Theorem [3l by linearity, we can 
assume / to be real-valued, and then, given r] > and (72 < / < fi'i as in the approximation 
theorem, we have 

|det(A„,)|E[(72(X„)] < |det(A„)|E[/(X„)] < | det(A„)|E[(7i(X„)], 

and hence 

< limsup I det(A„)|E[/(X„)] - liminf | det(A„)|E[/(X„)] < [ (g^ - g2)dm < r], 

n "'J 

which proves the result since r^ > is arbitrary. Similarly, the proof of (j5]) is performed in 
standard ways. D 
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Remark 3. To illustrate why our results are generalisations of the local theorems, let us 
analyse a particularly simple situation. We assume that d = 1 and that the random variables 
Xn have characteristic functions ipn such that fnit/bn) converge to ip{t) in -L-^(M), with 
bn — )■ +C)0 (such situations are related, but less general, than the classical results discussed 
in Section EH] or in [S] and [25].) In that case, the density functions /„ of X^/bn exist, are 
continuous and converge (in L^(M) and uniformly) to a continuous density function /. For 
a bounded interval (a,/3), we obtain 

6„P[X„ G (a, /?)] = 6„P[X„/6„ e {a/bn, /3/6„)] 

rl3/b„ 
= bn fn{x)dx ^ f{0){f3-a), 

Ja/bn 

by elementary calculus. 

It may be worth remarking explicitly that it is quite possible for this theorem to apply in a 
situation where the constant ^(0) is zero. In this case, the limit gives some information, but 
is not as precise as when the constant is non-zero. For instance, consider the characteristic 
function (p{t) = 1/(1 — it)"^, which corresponds to the sum Ei + E2 of two independent 
exponential random variables with density e~^dx on [0, +oo[. An easy computation shows 
that the density for (p itself is xe~^ (supported on [0, +oo[), and for X„ = n{Ei + E2), we 
have mod-if convergence with A„t = nt, leading to the limit 

hm nF[a < X„ < /3] = 

n— >+oo 

for all a < p. Note that any other limit c{/3 — a) would not make sense here, since Xn is 
always non- negative, whereas there is no constraint on the signs of a and /3... 

However, in similar cases, the following general fact will usually lead to more natural 
results: 

Proposition 1 {Mod-(f convergence and shift of the mean). Let d > 1 be an integer, and 
let (Xn) be a sequence of M.'^ -valued random variables such that there is mod—(f convergence 
with respect to the linear maps A„. Let a ^M.'^ be arbitrary, and let an G M'^ 6e a sequence 
of vectors such that 

(6) lim S„a„ = a, 

for instance an = Ana. Then the sequence Yn = Xn — an satisfies mod-ip convergence with 
parameters An for the characteristic function 

m = ^(t)e-^*-. 

In particular, for any continuous function f on M'^ with compact support, we have 

lim \detiAn)\E[fiXn-an)] = ^ia) f f{x)m{dx), 

where fi is the probability measure with characteristic function ip, and for any bounded 
Jordan-measurable subset B C W^ , we have 

(7) lim I det(A„)|P[X„ - a„ G 5] = ^[a)m{B). 

n^+00 dm 



Proof. This is entirely elementary: ip is of course integrable and since 

we have E[e*^'**'^"] = (y9„(S*t)e~**'^""", which converges locally uniformly to ip{t) by our 
assumption ([6]). Since the modulus of the characteristic function of Yn is the same, at any 
point, as that of X„, Property H3 holds for (Yn) exactly when it does for (X„), and hence 
mod-ip convergence holds. If / = dfi/dm, the density of the measure with characteristic 
function ip is g{x) = f{x + a), and therefore the last two limits hold by Theorem |H D 

In the situation described before the statement, taking a„ = en with c > leads to the 
(elementary) statement 

lim nF[a + en < Xn < l3 + en] = ee'^^S - a). 



hoo 



Even when the density of n does not vanish at 0, limits like ([7]) are of interest for all a 7^ 0. 
Another easy and natural extension of the local limit theorem involves situations where a 
further linear change of variable is performed: 

Proposition 2 (Local limit theorem after linear change of variable). Suppose that (X^) 
satisfies mod-Lp convergence relative to An and S„. Suppose that (Tn) is a sequence of linear 
isomorphisms T^: M'* -^ M'^ such that T~^ — > and S„T„ — ;■ 0. Suppose also that the 
following balancedness condition holds: there is a constant C such that |(S„T„)*t| < 1 implies 
that |S*t| < C. Then the sequence T^^Xn also satisfies the conditions of the theorem, and 
in particular for any hounded Jordan measurable set B we have 

l'^'="-^»)lplX„er„B]^^(0)m(B). 



|det(T„)| dm 

Proof. Let us put X„ = T~^Xn and (p{t) = E.[exp{it.T~^Xn)] = ipn{(T^^)*t). Clearly the 
sequence Sn^n tends to zero and S„T„(X„) = SnX„ tends to /i in law. The only remaining 
thing to verify is the uniform integrability condition. Let us look at 

(^„((S„,T„)*t)l|(s„T„)*t|<fc = V5„((S„,)*t)l|(s„r„)*t|<fc- 
Because of the balancedness condition we get that 

l|(S„T„)*t|<fc < l|(S„)*t|<CA:- 

The rest is obvious. D 

Remark 4. The balancedness condition is always satisfied if rf = 1. In dimension d > 1, there 
are counterexamples. In case the ratio of the largest singular value of S,„ to its smallest 
singular value is bounded, the balancedness condition is satisfied (this is an easy linear 
algebra exercise). See also [18] for the use of such conditions in mod-gaussian convergence. 
To see that for rf = 2 it is not necessarily satisfied take the following sequences: 



2.3. Conditions ensuring mod-ip convergence. We now derive other equivalent condi- 
tions, or sufficient ones, for mod-(/9 convergence. First of all, the conditions HI, H2, H3 have 
a probabilistic interpretation. We suppose d = 1 to keep the presentation simple. Instead of 
taking the indicator function l|E*t|<fc; we could have taken the triangular function A^ defined 
as Afc(O) = l,Afc(2A;) = = A(— 2A;), Ak{x) = for |a;| > 2k and A^ is piecewise linear 
between the said points. The function Ai is the characteristic function of a random variable 

Y (taken independent of the sequence X„). Hence we get that the sequence X„ satisfies 
HI, H2, H3 if and only if for each k > 1, the characteristic functions of Z„ = S„ (X„ + ^Y) 
converge in L^{M.'^) to (/?. Indeed the characteristic function of Z„ equals 9?„(S„t)Afc(E„t). 
There is no need to use the special form of the random variable Y. 

In fact we have the following: 

Theorem 5. Suppose that for the sequence X„ the conditions H1,H2 hold. The condition 
H3 holds as soon as there is a random variable, V , independent of the sequence X„ such 
that for each e > 0, ¥, [exp (itS„ (X„ + eV))] tends to cp in L^{W^). 

Proof. Let ■?/'(t) = E [exp {it ■ V)] be the characteristic function of V. The hypothesis of the 
theorem is equivalent to the property that for each e > 0, the sequence 

is uniformly integrable. Let 6 > he such that for |t| < 6, \'ip{t)\ > 1/2. Then the uniform 
integrability of the above mentioned sequence implies for each e > 0, the uniform integrability 
of the sequence 

\MKme\^-t^t\<5 < 2 \MKtMeKt)\ ■ 

This ends the proof. D 

Remark 5. We suppose that d = 1. For higher dimensions the discussion can be made along 
the same lines but it is much more tricky. Polya's theorem says that if 7: M+ — )■ IR+ is a 
convex function such that 7(0) = 1, lim.x^+oDl{x) = 0, then there is a random variable 

Y (which can be taken to be independent of the sequence X„, such that v^y(t) = 7(|it|). 
The characteristic function of S„(X„ + Y) is then 7(|S„t|)(/9„(S„t) and hence is a uniformly 
integrable sequence. Since ||Sn|| -^ 0, we see that T,n{Xn + Y) tends in law to /i with 
characteristic function ip. The convergence is much stronger than just weak convergence. In 
fact the random variables S„(X„ + y) have densities and because the characteristic functions 
tend in L^ to ip, the densities of S„(X„ + Y) converge to the density of /i in the topology of 
L\R). 

Remark 6. Adding a random variable Y can be seen as a regularisation (mollifier). Indeed 
adding an independent random variable leads to a convolution for the densities. In our 
context this means that the distribution of X„ is convoluted with an integrable kernel (the 
density of Y). The regularity of the law of Y is then passed to the law of X„ + Y. In prob- 
ability theory such a mollifier is nothing else than adding an independent random variable 
with suitable properties. 

We can go one step further and replace the condition for each A; by a condition where we 
use just one random variable. This is the topic of the next theorem. 

Theorem 6. Suppose that the hypotheses H1,H2 hold. Then H3 is also equivalent to either 
of the following: 



(1) There exists a non-increasing function 7: M+ — )■ R+, such that 7(0) = 1, < 7 < 1, 

lim 7(0;) = 

and such that the sequence 7(|S*t|)(y9(S*t) is uniformly integrahle. 

(2) There exists a non-increasing convex function 7: ]R_|_ — )■ M+, such that 7(0) = 1, 
<7 < 1, 

lim 7(3;) = 

x-^+oo 

and such that the sequence 7(|S*t|)(/9(S*t) is uniformly integrahle. 

Proof. It is quite clear that (1) or (2) imply H3, since for any A; > 0, we obtain 

|^„(s;t)|i|s;,i<fc < JL7(|s:t|)|^(s;t)|, 

and therefore the desired uniform integrability. 

For the reverse, it is enough to show that H3 implies (2), since (1) is obviously weaker. 
For a; > we define 

g{x) = sup ||v?„(S;t)l|s-t|<x+i||Li( 



The function g is clearly non-decreasing. Let us first observe that there is a constant 
C > such that for x big enough, exp (— f^ g{s) ds) < exp{—Cx). We define 



00 



7(0;) = a exp I — / g{s) ds j rf-u, 

where a is chosen so that 7(0) = 1 (since the integrals converge, this function is well defined). 
The function 7 is also convex and tends to zero at 00. Furthermore 

7(3;) < a —-— exp I — / g{s) ds j du < —— exp I — / g[s) ds 

Jx 9[^) \ Jo J 9[^> V Jo 

from which it follows that /g 'y{x)g{x) dx < 00. We now claim that 7(|S*t|)(p„(S*t) is 
uniformly integrable. Because the sequence is uniformly bounded we only need to show that 
for each e > there is a /c such that 

7(|s;t|)|v9„(s;t)|i|,|>,m(rft)<£ 

JR'' 

for all n > 1. 

For k, K integers, we split the integral as follows 



7(|S;t|)|v9„(S;t)|l|,|>,m(rft) < / 7(|S:t|)|v9„(S;t)|l|,|>a|s*i|<i.m(rft) 

+ / ^{\Kt\)Wn{Km\t\>kM^it\>Km{dt) 

+ / 7(|S:t|)|v9„(S:t)|l|s.,|>i,m(rft). 



The last term is dominated as follows: 

/ 7(|S:t|)|(^„(S:t)|l|s*t|>i.m(rft) < J2^{1) f \MKt)\midt) 

1>K 

,^^ \ Jo J 



1>K 



which can be made smaller than e/2 by taking K big enough. Once K fixed we use the 
uniform integrability of the sequence v9„(S*t)|l|E*t|<x and take k big enough so that we get 
for each n: 



|v9„(S;;t)|l|t|>fcl|2.t|<xm(rft) < e/2. 

This completes the proof. D 

In particular, we get a sufficient condition: 

Corollary 1. Suppose that the sequence Xn satisfies the following: 

(1) H1,H2 hold; 

(2) There is a non- decreasing function c: W ^ M_|_, c(0) = 1 as well as an integrable 
function h: R'^ -^ R+ such that \ipn{t)\ < h{A*j) c{\t\) for all n and t e R"^. 

Then the property H3 holds as well. 

Proof. This is clear from the previous theorem, since |v5n(^)| ^ h{A^t) c{\t\) for all t implies 
that for all t we have 

l^"(^n^)l < hit) 

which verifies (1) in Theorem O □ 

3. Applications 

In this section, we collect some examples of mod-ip convergence, for various types of limits 
if, and therefore derive local limit theorems. Some of these results are already known, and 
some are new. It is quite interesting to see all of them handled using the relatively elemen- 
tary framework of the previous section. The coming subsections are mostly independent of 
each other; the first few are of probabilistic nature, while the last ones involve arithmetic 
considerations. 

3.1. The Central Limit Theorem and convergence to stable laws. In this section we 
suppose that (X„)„>i is a sequence of independent identically distributed random variables. 
The central limit theorem deals with convergence in law of expressions of the form ^i+-+-^" _ 
a„, where 6„ are normalising constants. We will suppose without further notice that the 
random variables are symmetric so that we can suppose a„ = 0. The possible limit laws 
have characteristic functions of the form exp(— c|up), where < p < 2 and where c > 0. 
For information regarding this convergence we refer to Loeve |2D]. The basis for the theory 
is Karamata's theory of regular variation. In this section we are interested in expressions of 
the form lim6„P[(Xi + . . . + X„) G B] for suitably bounded Borel sets B. 

10 



For the case E[X^] < cxo, the problem was solved by Shepp [23]. The multidimensional 
square integrable case was solved by Borovkov and Mogulskii [1] and Mogulskii [22]. The 
case p <2 was solved by Stone ^25j and at the same by BretagnoUe and Dacunha-Castelle [5] 
(see also Ibragimov and Linnik [I2|)- Such theorems are known as local limit theorems. 

Theorem 7. Suppose that the non-lattice random variable X is symmetric and that it is 
in the domain of attraction of a stable law with exponent p. More precisely we suppose that 
Xi+...x„ converges in law to a probability distribution with characteristic function exp(— |t|P), 
< j9 < 2. Then for Jordan-measurable bounded Borel sets, we have 

lim 6„P[(Xi + ... + Xn)eB] = Cp m{B), 

where Cp = ^ J_^ exp(— |t|^) dt. Suppose moreover that < r„ — )■ +oo in such a way that 
^ -)■ +CX), then 

lim ^P \ Xi + ...+Xn ^ 1 ^ Cpm{B). 

n-5-+oo Tn T„ 



In order to prove this theorem, we first observe that, when HI and H2 are satisfied, the 
condition H3 of uniform integrability is equivalent with classical conditions that arise in the 
current context. 

Theorem 8. Under the hypotheses H1,H2, the hypothesis H3 is equivalent to the validity 
of the following two conditions: 

• H3'. For all k > e > , we have 



lim |det(A„)| / |^„(t)| m(rft) = 0. 

^^+°° Je<\t\<k 



H4'. For all r] > 0. there is a > 0,6 > such that 



limsup / |v9„(S*t)| m((it) <t]. 

n-^+oo J a<\t\:\T.tt\<e 



n-^+oo Ja<|t|;|S*i|<£ 

Proof. First suppose that H3 holds, i.e., for each A; > 0, v9n(S*t)l|E*i|<fc is uniformly inte- 
grable. Since E„ — )■ 0, we immediately get 

lim/ \det{An)(pnit)\m{dt) = |v9„(S*s)| m(rfs) -^ 0, 

" A Je<\T.^s\<k 

which is H3'. To establish H4', let us first remark that (using H2) we have 

lim(^„(S;t)l|E.t|<fc = v{t) 
for all t. Then we take a > such that for given 77 > we have 

\(p{t)\m{dt) <r]. 

\t\>a 

Take now e > and observe that by uniform integrability 

lim/ |(y9„(S*t)| m((it) = / \ip{s)\7n{ds) < 7]. 

"■ "'a<|f|;|E*t|<£ -''|s|>a 
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Now we proceed to the converse and we suppose that H1,H2,H3',H4' hold. We first 
show that given 77 > 0, the sequence v9n(S*t)l|s*t|<A; has up to 77 all its mass on a ball of 
radius a. Given r^ > we can find a, e > such that 



lim / |v?„(S*t)|m((it) < 77. 

" J a<\t\:\T.lt\<e 



'«<|t|;|S*t|<e 

Then according to H3' and H4', we can find hq such that for all rz > rzo we have 

V9„(S;t)|m(rft)<2r/, 

a<\tl\J]'„t\<e 



\v{Kt)\m{dt)<r^. 

'e<|S*t|<fc 

Increasing a allows us to suppose that the same inequalities hold for all ri > 1. So we get 

that 



|^(S:t)|m(rft)<3r/. 

'|S*t|<fc;|i|>a 

Since the sequence is uniformly bounded we have proved uniform integrability. D 

Proof of Theorem [^ We have here </?„ = ip^ where tp is the characteristic function of a 
random variable in the domain of attraction of a stable law. Property H3' follows since the 
sequence ipn tends to zero exponentially fast, uniformly on compact sets of ]R^\{0}. Moreover, 
Property H4' is known as Gnedenko's condition (see Gnedenko and Kolmogorov [8] or the 
discussion of I2 (resp. J3) in Ibragimov and Linnik [121 p. 123]) (resp. [T?', p. 127]). Thus 
the hypotheses HI, H2, H3', H4' are fulfilled in this setting. D 

Remark 7. (1) The proof of Property H4' is based on the regular variation of ip around 0. 
The fact that regular variation is needed suggest that it is difficult to get a more abstract 
version of this property. 

(2) Taking the most classical case where p = 2 and (X„) independent and identically 
distributed, it is easy to check that the stronger condition H2' (i.e., (jl]) is not valid, except 
if the Xn are themselves gaussian random variables. Thus the setting in this paper is a 
genuine generalization of the original mod-gaussian convergence discussed in p^ . 

3.2. The winding number of complex Brownian motion. We take a complex Brownian 

Motion W , starting at 1. Of course we can also see W a.s a, two-dimensional real BM. The 
process W will never attain the value and hence, by continuous extension or lifting, we can 
define the argument 6. We get Wu = RuQy^'l^iiOu) where 6*0 = and i?„ = \Wu\- The process 
6 is called the winding number, see [2T] . Spitzer in [21] computed the law of 6^ and gave its 
Fourier transform, and a more precise convergence result was given in [3]. 
The characteristic function is given by 

E [exp(zt^J] = i^f" (i-) exp (-i-) (/(„_,)/, (i-) + /(|,^,)/, (i- 



where Iu{z) denotes the /-Bessel function, which can be defined by its Taylor expansion 

1 /Z\''+2m 

.2, 






m>0 
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Using elementary properties of Bessel functions, Spitzer deduced that ^^ converges to a 
Cauchy law with characteristic function ip{t) = exp(— |t|) and density ^jr^- 

Theorem 9 (Mod-Cauchy convergence of the winding number). For any sequence (u„) of 
positive real numbers tending to infinity, the sequence X„ = Ou^ satisfies mod-ip convergence 
with d=l, <f{t) = exp(-|t|), An{t) = Al{t) = (\ogUn)t/2. 
In particular, for any real numbers a < b, we have 

lim i^^ F[a<eu<b] = -{b-a). 

li— >oo 2 TT 

Although this is a very natural statement, we have not found this local limit theorem in 
the literature. 

Proof. The conditions HI and H2 of mod-y? convergence are clear, the second by Spitzer's 
Theorem. To check the uniform integrability condition H3, we take k > 1 and we proceed 
to bound 

for t > 0. But if |S;t| < k, we have 

1 ^ |s;t|±i ^ k + i 

2 - 2 - 2 
and < 7^ < 1 for n large enough. The Taylor series expansion shows immediately that 
there exists Ck > such that 

\Iu{z)\<Ck 
uniformly for u real with —1/2 < v < -^ii and z E C with \z\ < 1, so that for |S*t| = 
2|t|/(logM„) < k, we have 

\iPn{Kt)\ < Bku;,'/' < Sfcexp(-^) 
where B^ = Cfc(7r/2)^/^. This gives the desired uniform integrability, in the form ([3]). D 

3.3. "Relaxed" Poisson variables. We present here a special case of a phenomenon which 
is related to Poisson approximation and therefore probably quite general: if P„, n > 1, de- 
notes a Poisson-distributed random variable with parameter An going to infinity, the sequence 

P _ \ 



X.r, 



,1/3 

An 



satisfies mod— (y9 convergence with (i = 1, ip{t) = e ' (i.e., for a standard gaussian) and 
Ant = Xn t. Indeed, H2 holds because 

ipnit) = e-*""''*exp(An(e**/""'' - 1)) = exp(-Ay=^tV2) exp(-2tV6)(l + o(l)) 

. ,, 1/3 

as n tends to infinity. Moreover, the next term in the expansion of the exponential e ' " 
shows that 



^"1;^ 



An. 



and the uniform integrability condition H3 therefore holds even for the range |t| < An . 
(Except for H3, this example was considered in [T71 Prop. 2.4].) 
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As a consequence, we get the local limit 

lim A^SpiaAy^ + A„ < P„ < hX"/' + A„] 



^2tt 
for any fixed a < b. 

Remark 8. Using the formula [T71 (4.9)], we see that the same mod-ip convergence property 
holds when P„ is replaced with a;„ defined as the number of cycles in the decomposition in 
cycles of a uniformly chosen random permutation in the symmetric group on n letter, with 
A„ = logra. These are well-known (see, e.g., [H [T7]) to be well-approximated by Poisson 
variables with these parameters. 

3.4. Dedekind Sums. In this section we give an application to Dedekind sums. Our limit 
theorems are based on the estimates in Vardi's paper [2S]- Let us recall the definition of 
Dedekind sums. We recall the standard notation 

[a;J = supjrT, G Z | n < x}, 

((x)) =x-[x\ -1/2, if x^Z 

= if X G Z. 

For natural numbers < d < c with gcd{d, c) = 1, the Dedekind sum is defined as 



^<-'^s,.eo 



0<A:<c 

For every A^ G N we define the finite probability space: 

^N = {{d, c) \ < d < c < N; gcd(rf, c) = 1}, 

#^ 
Pat [A] = the normalised counting measure, 

X7v((i, c) = s{d, c). 
The distribution of Xn is symmetric as easily seen by using the measure preserving trans- 
formation {d, c) ^ {c — d, c). It is well known that jj^Vt^ / i ^^ ) — )■ 1, see e.g section 3.4 of 
this paper. Vardi [261 Prop. 2] proved an asymptotic formula which implies the following: 
Proposition 3. For < |t| < 1/4 we have that: 

|<^Ar(27rt)| <CiV-l*l+o(iV-i/=^), 
where C is an absolute constant and where the last term is uniform in t. 

Remark 9. The result of [26] actually gives the same result for larger values of t, but the 
error term is only smaller than the main term when |t| < 2/3. 

As a consequence of the same proposition in [26], we get that for t ^M.: 

^^ ( bi^ ) ^ ^^^^ ^ exp(-|t|), 

the characteristic function of a standard Cauchy random variable with density ^i-^^2\ - 

The bound given by Vardi does not allow to show a mod—*/? (in this case "mod-Cauchy") 
convergence, but it suffices to obtain the following weaker statement: 
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Proposition 4. For any sequence (rjy) such that t^ — )■ +00 and — )■ +00, the 



TN 



t. 



satisfies mod-ip convergence with Ajvt 
Hence, for every bounded Jordan-measurable set B gM., we have 



logN 
27rrjv 



TJV 



sequence 



logN 

271TN 



F 



N 



X 



N 



tn 



eB 



-)■ -m{B). 

71 



Proof. We only have to show that for each k, the sequence 

(PNi2nt/ log N)l 2.\t\rr, ^j^ 

log N — 

is uniformly integrable. This is seen as follows: if ^^ J^ < k and if A^ is big enough, then 
|t|/logA^ < 1/2. Consequently for A^ large enough (depending on k), we get 



(PN (27rt/ log A^) 1 2n\t\Tj^ 



logiV 



<k 



<Cexp(-|t|)+V^^(t), 



D 



where ||V'Af||i ^ CN ^^^logN. This implies uniform integrability of the sequence. 

Remark 10. The more precise local limit theorem 

logiV 1 

Zn TT 

is in fact valid, as proved by Bruggeman [6]. Our methods do not seem to lead to this result 
only using mod-cp convergence. 

3.5. The (^-distribution. The (^—distributions are purely atomic, infinitely divisible, prob- 
ability distributions, denoted fi'^, which were considered by Khintchine and studied in more 
detail in \19\. 

The measure /x'^ is defined for a > 1 as the measure supported on the points {— log(r2) | 
ra G N; ra > 1}, such that 



/i" 



log n) 



n 



C(^) 



for n > 1. Its characteristic function is then given by 



<^ 



^(^) = E 






n>l ^^ ' 

The limit of interest here is when cr J, 1. Since the zeta function can be written 



(is) 



s-V 



where C*{^) defines an entire function of s G C (i.e., the zeta function has only a simple pole 
with residue 1 at s = 1), the behavior of 'p'^{t) is easy to understand, namely 

Cia + tt) 1 C{c^ + ^t) 



^"{t) 



CM 



1 + ^ 



CM 



Thus, if X'^ are random variables with law /i'^, we see that (cr — 1)X" converges in law to a 
'negative" exponential distribution supported on ] — 00, 0] with density e^. The characteristic 
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is not integrable, hence we can not apply our results. To work around this, we consider 
independent copies X{^,X^ of random variables having the law /i°", and define 

F<^ = XI - X^. 
These random variables have characteristic function given by 

|C(a + ^t)P 1 \C{^ + ^t)? 



w^m 



a^y i + (^ ^*(^) 



and hence {a — l)y°" converges in law, as cr J, !> to a double exponential (or Laplace) dis- 
tribution, with characteristic function ip(t) = j;^^ and density |exp(— |x|). Thus conditions 
HI and H2 are now satisfied (in the version of a continuous limit a J, 1). Moreover, if 
1 < (T < 2 and t ranges over the set where \{a — l)t\ < k, for A; > fixed, the values of 

\C{a + ^t{a - 1))\' 



C(a)2 

vary in a bounded set. This shows that H3 also holds, and we can apply Theorem |H it 
follows that 

lim F[a <¥"" <b] = -{b- a). 

0-4.1 a — 1 2 

for all —00 < a < b < +00. We can make this limit explicit: indeed, V^ takes values of the 
form \og{k) — log(?2) = log{k/n) where k,n > 1. The probability that V^ = \og{k/n) for 
{k, n) = 1 (i.e. k and n are coprime) is easily seen to be 



„>. «")= «-) 



Hence the limit becomes, for < a < /3, the formula 

C(2a) 



E '---A'^^i^ 



{k,n)=l 
a<k/n<P 

as cr 4, 1, which is equivalent (since C(2) = 7r^/6 and C,{a) ~ (a — 1)) to 

(.-1) 5: r",r"^-ilog(g. 

(fc,n)=l ^ ^ 

a<k/n<li 

We could not find any reference to this statement, so it might be new (although it could 
certainly be proved with more traditional methods.) 

3.6. Squarefree integers. This section is motivated by a recent paper of Cellarosi and 

Sinai [7J, who discuss a natural probabilistic model of random squarefree integers. As we 
will see, some of its properties fall into the framework of mod-(^ convergence, with a very 
non-standard characteristic function ip. 

The set-up, in a slightly different notation than the one used in [7], is the following. We 
fix a probability space Vt that is big enough to carry independent copies of random variables 
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rip, with index p running over the prime numbers, with the following distribution laws: 
F[np = 1] = F[rip = -1] = -^—, nVp = 0] ^ ^ ^ 



We consider the random variables 

^n = ^ Vp^ogp, Qn = exp(X„) 

'P<Pn 

for n > 1, where Pn is the n-th prime number. 

The link with [7] is the following: in the notation of [T, Th. 1.1], the distribution of X„ is 
the same as that of the difference 

(Cn-C)logPn 

of two independent copies Cn and ('^ of the random variables variables 

P<Pn 

of [3 Th. 1.1], where the Up are independent Bernoulli variables with 

n^p = 0] = -^, n^p = 1] ^ 



p+i ' " ' p+i 

These random variables Up are very natural in studying squarefree numbers. Indeed, a 
simple computation shows that Up is the limit in law, as x — )■ +oo, of the Bernoulli random 
variables z/p .j. defined by 

|. \{n < X \ X squarefree and divisible by p}\ 

^'^ \{n < X \ X squarefree}! 

for fixed p. 

By definition, the support of the values of exp{(n) is the set of squarefree integers only 
divisible by primes p < Pn, and for Qn, it is the set of rational numbers x = a/b where a, 
b > 1 are coprime integers, both squarefree, and both divisible only by primes p < pn- It 
is natural to see them as giving probabilistic models of these numbers. We obtain mod-y? 
convergence for X„: 



Theorem 10. Let 



^M = exp(-4^'sm=(|)^) 



for t G M. Then ip is an integrable characteristic function of a probability distribution on M, 
and the sequence (X„) satisfies mod-ip convergence with d = \ and Anit) = v4* (t) = (logp„)t. 

The proof is quite similar in principle to arguments in [3|, though our presentation is more 
in the usual style of analytic number theory. 
We start with the easiest part of this statement: 

Lemma 1. We have ip G L^(]R), and in fact 

(8) Ht)\<C\tr 

for \t\ > 1 and some constant C > 0. 
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Remark 11. The characteristic function of the hniit in law of the (non-symmetrised) random 
variables (n used in [T^ Th. 1.1] only decays as t~^ when |t| — )■ +00, and hence is not 
integrable, which prevents us from applying our results directly to those variables. As we 
will see, this is quite delicate: changing the constant 4 to a constant < 2 would lead to a 
failure of this property. 

Below, we will see that Theorem H] is not valid for the variables (logj9„)Cn- 

Proof. Integration by parts gives that 

/"* sin^ X 1 

/ dx = -\ogt + b(t) 

Jo X 2 

where b{t) tends to a constant for t — t- 00. From here we deduce that 



4/ sin2('— ")— = 21og|t| +c(t) 
Jo \ 2 / V 



where c(t) remains bounded. As a result we get ([8]), which proves that (p E L} since the 
function is continuous. (Alternatively, one can check that 

(9) ¥?(t) =exp(-27-21og|t| +2Ci(t)) 

where 7 is the Euler constant and Ci(t) is the cosine integral function, and use the properties 
of the latter.) D 



Proof of TheoremUU Let 



Y - ^ X 

logPn 



and let ipn be the characteristic function of y„, which we proceed to compute. 
With X = a;„ = p„, we have first 



^„(t) = E[exp(^ty;)] = JJe 



logp 
logx 



exp ( ^ i__ ^ trjp 



p<x 

nf p'^ + l 2p / logp 

TT fi - ^^ f 1 - cos f tj^ 

nfi-^^sin^f^j^ 



(10) =expK^log(l 



4p . 2 ftlogp 
sm - 



for all t G M. Now we assume t 7^ (since for t = 0, the values are always 1). We first show 
pointwise, locally uniform, convergence. 
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The idea to see the hmit emerge in the sum over p is quite simple. First of all, we can 
expand the logarithm in Taylor series. We have 

2/tlogp' 



lim 7 > p 



k>2 p<x 



sm 



21ogx 



0, 



for t in a bounded set, by dominated convergence. This allows us to restrict our attention 
to 

tlogp' 



11 



4^p Sin^[ 



p<x 



21ogx 



(we also used the fact that 4:p/{p + 1)^ is equal to A/p up to terms of order p~^.) Now, for 
P ^ y, where y < a;^'''*' is a further parameter (assuming, as we can, that this is > 2), we 
have also 



\^P ^sin2(^ 



p<y 



tlogp 
21ogx 



< 



21ogx 



J2p '(logp)'<t' 



p<y 



(logx)^ 



Thus, for if we select y = y{x) < x^''*' tending to infinity slowly enough that logy 
o(logx), this also converges to as x — )■ +oo, and what remains is 

. V^ -1 • 2/^logP 
-4 2_. P sm 



y{x)<p<x 



21ogx 



We can now perform "back-and-forth" summation by parts using the Prime Number The- 
orem to see that this is 



-4 



-1 . 2i"t^ogu\ du 
u sm [ —, jj— h o(l) 



ly{x) v21ogX/logM 

as X —7- +00 (apply Lemma [2] below with B = 2 and with the function 

tlog-u 



f{u) = -sin 

u V 2 log X 



with 



fin) 



sm 



which satisfies 
(12) 



t\ogu\ t_ 

w" \2\ogxJ M^(logx) V21ogx 

\f{u)\<u-\ \f'{u)\<(l + -^) 



t log u\ ft log u 

sm I — ) cos — 

V 2 log X 



u 



the integral error term in Lemma |2] is then dominated by the tail beyond y{x) of the con- 
vergent integral 

^ du 



12 M(logM)2' 

and the result follows) . Performing the change of variable 

logM 

v = -, , 

logx 
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we get the integral 

, f^ ■ 'i/'tv\dv 

-4 / sm^ — — , 

Jlog{yix))/ilogx) \2 J V 

which converges to (f{t) as x — )■ +00. 

To conclude the proof of Theorem [TOl we will prove the following inequality, which guar- 
antees the uniform integrability condition H3: for any k > 1 and t, n with |t| < A;(logx) = 
k(\ogpn), we have 

(13) Mt) < |<^(t)|exp(Cloglog3|t|) 

which gives the desired result since we know from (|H]) that ip decays like |t|~^ at infinity. 

We can assume that |t| > 2. Now we start with the expression fIlOp again and proceed to 
deal with the sum over p < x in the exponential using roughly the same steps as before. To 
begin with, we may again estimate the sum f lTT|) only, since the contribution of the others 
terms is bounded uniformly in t and x: 

k>2 p<x *= fc>2 p 

which is a convergent series. After exponentiation, these terms lead to a fixed multiplicative 
factor, which is fine for our target (IT3|) . 

We next deal with the small primes in flTTj) : since |t| < klogx, the sine term may not lead 
to any decay, but we still can bound trivially 

E_-, . r,/t\ogp\\ -^-^ _-. , , 
P sm — < Z^P <loglog?/ 
V / log X / I -^^ 
p<y p<y 

for any y < x (by a standard estimate). We select y = \t\ > 2, and this becomes a factor of 
the type 

exp (Clog log t) 

(after exponentiating), which is consistent with (TT3l) . 

We now apply Lemma [2] again, writing more carefully the resulting estimate, namely 

A \^ ~1 ■ 2ft^Ogp\ r _i . 2/tl0gM\ du , ^f l + k 

— 4 > p sm — = — 4 / u sm — h O — — 

^—^ V21ogx/ Jy V21ogx/log'U \[\ogy)^ 



iiogv/iogx ^2/ V(logy)^ 



(using the bound flT^ ). with an absolute implied constant. The remainder here is again fine, 
since y = \t\ > 2 by assumption. 

Now, to conclude, we need only estimate the missing part of the target integral (which 
runs from to 1) in this expression, namely 



i 



logS//logx _ ^^ 

V sin [ — \dv. 



'0 

We write 

flog y/ log x /tV\ /"l*l"^ r^ogy/logx 
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where the first terms is bounded by 



r\t\-^ 
{t/2f I vdv < 1, 
Jo 



and the second by 

flog y/ log X 



Mogy/iogx ^ f\t\\ogy\ 

/ v~^dv = \og[— ) < logfA; logy) = logfA; log Itl). 

i|t|-i ^ logx / 

Putting the inequalities together, we have proved (fT3|) . and hence Theorem [TOl D 

Here is the standard lemma from prime number theory that we used above, which expresses 
the fact that for primes sufficiently large, the heuristic - due to Gauss - that primes behave 
like positive numbers with the measure du/{\.ogu) can be applied confidently in many cases. 

Lemma 2. Let y > 2 and let f be a smooth function defined on [y,+oo[. Then for any 
A > 1, we have 

^^^ Jy logM V(logx)^ (logy)^ Jy (logw) 

where the sum is over primes and the implied constant depends only on A. 

We give the proof for completeness. 

Proof. We use summation by parts and the Prime Number Theorem, which is the case 
f{x) = 1, in the strong form 

f^ du / X 

T^{x) = I V0{- -J 

J2 logu V(log2;)^ 

for X >2 and any A> 1, with an implied constant depending only on A (this is a consequence 
of the error term in the Prime Number Theorem due to de la Vallee Poussin, see e.g. [HI 
Cor. 5.29]); this leads to 






{u)'K{u)du, 

y<p<si 



and after inserting the above asymptotic formula for 7r(a;) and vr('u), we can revert the 
integration by parts to recover the main term, while the error terms lead to the result. D 

We now derive arithmetic consequences of Theorem [TUJ Applying Theorem 2, we get 

(14) lim (logp„)P [Xn e]a, b[] = {h - a)r] 

n— >'+oo 

where 

r^ = - r^exp (-A f'^^dv] dt = - /^"^ e^^^'^"^)^, 
' 2nJ_^ "^y J, V J 2nJ_^ |t|2' 

the last expression coming from (|9]). Using the relation between (f and the Dickman-de 
Bruijn function p, namely 

v?(t) = mn-t) 
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where ilj{t) is the Fourier transform of e "'p{u) (this follows from [7| Th. 1.1, p. 5]), one gets 



1] = e '^^ 



I p{ufdu = 0A5A867 . 



^n — - 



(the numerical computation was done using Sage). 

This arithmetic application could certainly be proved with more traditional methods of 
analytic number theory, when expressed concretely as giving the asymptotic behavior as 

n —7- +00 of 

E f 

a<r/s</3 

but it is nevertheless a good illustration of the general probabilistic framework of mod-y^ 
convergence with an unusual characteristic function. 

Although our theorem does not apply for the random model of [7] itself, it is quite easy to 
understand the behavior of the corresponding probabilities in that case. Indeed, denoting 

Yn = exp((logp„)C„), 
which takes squarefree values, we have 

^ fc<e" 

for any fixed a G M, where a; = p„, /i^(/c) is the indicator function of squarefree integers and 
Zx is the normalizing factor given by 

For X large enough and a fixed, the second condition is vacuous, and hence this is 

As observed in [7| (3)], we have Z^ ~ e'^({2)~^ logx, and hence we get 

hm (logp„)P[r„ < el = C(2)e-^ V ^. 

fc<e" 

When a is large, this is equivalent to e~'^a (another easy fact of analytic number theory), 
which corresponds to the local limit theorem like (1141) . but we see that for fixed a, there is 
a discrepancy. 

There is one last interesting feature of this model: the analogue of Theorem [10] for polyno- 
mials over finite fields does not hold, despite the many similarities that exist between integers 
and such polynomials (see, e.g., [I7j for instances of these similarities in related probabilistic 
contexts.) 

Precisely, let g > 1 be a power of a prime number and F^ a finite field with q elements. For 
irreducible monic polynomials vr G IFq[X], we suppose given independent random variables 
rJT^, T]'^ such that by 

= ±1] = m = ±1] = (^^, nv. = 0] = m = 0] = ^^^ 
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where |7r| = q^'^si->r\ Then for ra > 1, let X„ be the random variable 

dcg(7r)<n 

where the sum runs over all irreducible monic polynomials of degree at most n. Then we 
claim that HI, H2 hold for X„, with the same characteristic function ip{t) as in Theorem [TOl 
and Ant = nt, but there is no mod— y? convergence. 

This last part at least is immediate: H3 fails by contraposition because the local limit 
theorem for 

lim nF[a < X„ < 6] 

n— >+oo 

is not valid! Indeed, X„ is now real-valued, and if ]a, 6[nZ = 0, the probability above is 
always 0, whereas the expected limit [b — a)r] is not. 

We now check H2 in this case. Arguing as in the beginning of the proof of Theorem [TUl 
we get 

dcg(7r)<n ^ VI I I y \ / / 

Expanding the logarithm once more, we see that it is enough to prove that (locally uni- 
formly in t) we have 



lim exp I —4 > - — - sin^ I 

\ deg(7r)<n 



^'^(")'^> ' - ^(t). 



2n 
We arrange the sum according to the degree of vr, obtaining 



where T^q{j) is the number of monic irreducible polynomials of degree j in Fg[X]. The 
well-known elementary formula of Gauss and Dedekind for T\.q{j) shows that 

n,(j) = ^ + 0{q^l^) 
for q fixed and j > 1, and hence we can write the sum as 

v^ 1 . 2 /deg(7r)t\ ^1 . 2 f 3t\ , n(^ -i/2 • 2(^3 



^___ I ■ - 1 / g ■" sm 

TT \ 2n J ^-^ i \2n I \^-^ \2n 

deg(7r)<n ' ' ^ ^ i=l "^ ^ ^ \i=l ^ 

As n goes to infinity, the second term converges to by the dominated convergence 
theorem, while the first is a Riemann sum (with steps l/n) for the integral 

. /'tv\ dv 
sm^ — — , 

and hence we obtain the desired limit. (This is somewhat similar to [Ij Prop. 4.6].) 
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Remark 12. A more purely probabilistic example of the same phenomenon arises as follows: 
define 

n 

i=i 
where {Dj,Ej) are globally independent random variables with distribution 

nEj = j] = nD, = j] = J, nE, = o] = p[d, = o] = i - ^. 

Then the sequence (X„) also satisfies HI and H2 for the same characteristic function ip{t) 
(by very similar arguments), and does not satisfy H3 since X„ is integral- valued. 

3.7. Random Matrices. Some of the first examples of mod-Gaussian convergence are re- 
lated to the "ensembles" of random matrices corresponding to families of compact Lie groups, 
as follows from the work of Keating and Snaith [T3], |15]. Using this, and our main result, 
we can deduce quickly some local limit theorems for values of the characteristic polynomials 
of such random matrices. 

We consider the three standard families of compact matrix groups, which we will denote 
generically by G„, where G is either U (unitary matrices of size n), USp (symplectic matrices 
of size 2n) or SO (orthogonal matrices of determinant 1 and sizqj 2n). In each case, we 
consider G„ as a probability space by putting the Haar measure /i„ on G„, normalized so 
that /i„(G„) = 1. The relevant random variables (X„) are defined as suitably centered values 
of the characteristic polynomial det(T — Qn) where Qn is a G„- valued random variable which 
is /i„-distributed. Precisely, define 



ar. 



if G = [/, 

I log(7™/2) if G = USp, 
I log(87r/n) if G = SO, 



and consider X„ = logdet(l — gn) — ctn', this is real-valued except for G = f/, in which case 
the determination of the logarithm is obtained from the standard Taylor series at z = 1. 
Now define the linear maps 



An{t) 



(log(l)) t otherwise. 



and their inverses E„ (these are diagonal so A* = An, S* = S„). 

Finally, let ip be the characteristic function of a standard complex (if G = U) or real 
gaussian random variable (if G = USp or SO); in particular HI is true. It follows from 
the work of Keating and Snaith that in each case (/?„(E„t) converges continuously to (p(t), 
i.e., that H2 holds. In fact, in each case, there is a continuous (in fact, analytic) limiting 
function $g(^) such that 

y.„(t) = y.(A;t)$G(t)(i + o(i)) 



The odd case could be treated similarly. 
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for any fixed t, as n goes to infinity. These are given by 



if G = f/, 



G{l+iti) 
G(l/2+it) II *cr - DCy, 

in terms of the Barnes G-function. Detailed proofs can be found in [HI §3, Prop. 12, Prop. 
15], and from the latter arguments, one obtains uniform estimates 

for all t such that |t| < n^/^, where C is an absolute constant. This immediately gives the 
uniform integrability for v5„(S*t)l|s„t|<fc since |S5^t| is only of logarithmic size with respect 
to n. In other words, we have checked H3, and hence there is mod-y^ convergence. 

Consequently, applying Theorem HI we derive the local limit theorems (already found 
in [H]): 

Theorem 11. For G = U, USp or SO, for any bounded Jordan-measurable set i? C M or 
C, the latter only for G = U, we have 

m( R^ 

lim I det(A„)| fi^ig e G„ | logdet(l - g) - an e B) = 

n-^+oo (ZTTj ' 

with d = 2 for G = U and d = 1 otherwise. 

Theorem [21 stated in the introduction, is the special case G = U, enhanced by applying 
Proposition [H 

As in [18, §4], one can derive arithmetic consequences of these local limit theorems, involv- 
ing families of L-functions over finite fields, by appealing to the work of Katz and Sarnak. 
The interested readers should have no difficulty checking this using the detailed results and 
references in |18j . 

Instead, we discuss briefly a rather more exotic type of random matrices, motivated by 
the recent results in |T6] concerning certain averages of L-functions of Siegel modular forms. 
In [ISl Rem. 1.3], the following model is suggested: let G„ = S'02„(M), with Haar measure 
fin, and consider the measure 

^nig) = - det(l - g)dfinig) 

on Gn- The density det(l — g) is non-negative on (?„ (because eigenvalues of a matrix 
in S0{2n,R) come in pairs e*^ e''^ , and (1 - e'^){l - e'^% > 0); the fact that this is a 
probability measure will be explained below. In probabilistic terms, this is the "size-biased" 
version of /i„. 

Theorem 12. Let X„ = logdet(l — gn) — |log(327rn), where gn is a Gn-valued random 
variable distributed according to Vn- Let (p be the characteristic function of a standard real 
gaussian. Then we have mod-ip convergence with Ant = (log^)^/^t, and in particular 



lim ./log ^ UnigeGnl logdet(l -g)-l log(327m) G B^ ""^^^ 



n^+oo v 2 y 2 / v27r 
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Proof. The characteristic function of ¥„, = log det(l — ^„) is half of the value at s = 1 + it of 
the Laplace transform E[e*'°§'^''*^^~^"''], where gn is Haar-distributed. The latter is computed 
for all complex s in [HI (56)], and we get 

ororJin^l _ jn'!jl+it)logdct{l-gn)-\ _ r,2n{l+it) TT ^(j + TO - l)r(j + jt + 1/2) 

^ ^" ^ ^" ii r(j-i/2)r(j + ^t + n) ■ 

At this point, the reader may check easily (by recurrence on n if needed) that this gives 
the right values E[e**^"] = 1 for t = 0, confirming the normalizing factor 1/2 used in the 
definition of Un- 

To go further, we transform the right-hand side into values of the Barnes function G{z), 
as in [m §4.3], to get 

n^Utv.^ = n2n(i+it) Gil/2) Gi2n)G{n + 3/2 + tt)Gil + tt + n) 
^ ^ G{3/2 + it) G{n)G{n + l/2)G{2n + l + it) ' 

Applying r{z)G{z) = G{z + 1), we transform this into 

2EUt^ni ^ o2n(wt) G{l/2) r(zt + n)r(zt + n + l/2) 
^ ^ G{3/2 + it) T{it + 2n) 

G{2n)G{n + it)G{l/2 + it + n) 



G{n)G{n + l/2)G{2n + it) 



and the last ratio of Barnes functions (together with the factor 2^™*) is exactly the one 
handled in jTSl Prop. 17, (4)]. With the asymptotic formula that follows, the Legendre 
duplication formula and r(l/2) = y/n, we deduce 

^, ,,y, G(3/2) T(2it + 2n) fn\-tV^ fSirY^^ ,^ 
2E[e'*^"] = 2 ) ' . , T,V r(-) { — ] (l + o(l)) 



G{3/2 + it) T{it + 2n) \2J \nj 

uniformly for |t| < n^^^. Since 

r{2it + 2n) _^^.^, 



r{it + 2n) 
in this range, we get 



(2n)^*(l + o(l)) 



E[e^*^1 = G||T^^(^"^)^^2vrn)^*/2(l + o(l)), 

uniformly for |t| < n^/^, and the result follows. D 

The most obvious feature of this exotic model of orthogonal matrices is the "shift" of the 
average; whereas, for Haar-distributed g G S'02„(M), the value logdet(l — g) is typically 
small (mean about log a/Stt/w) , it becomes typically large (mean log v327m, of similar order 
of magnitude as the mean for a symplectic matrix of the same size) when g is considered 
to be distributed according to Un- This is consistent with the discussion in [121 Rem. 1.3], 
especially since the "limiting function" that appears here is ^usp- 
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3.8. Stochastic model of the Riemann zeta function. The following "naive" model of 
the Riemann zeta function on the critical line is surprisingly helpful. The basic ingredient 
is a sequence of iid variables 1^ : f2 — )■ T where T is the unit circle in C and the variables 
Yp are uniformly distributed over T. For notational ease the sequence is ordered by the 
prime numbers. In what follows p will always denote a prime number. The random variables 

we consider are constructed as follows. First we take finite products Z„ = Y[p<n [^ — rj- 

If we replace the factors Yp by exp{ipt), then the product appears in the study of the 
Riemann (^-function. An easy application of Weyl's lemma on uniform distributions shows 
that (exp(ipt))^^ defined on [0,T] (with normafised Lebesgue measure) tend (as T — )■ oo) 
to {Yp)p<n- The random variables X„ are then defined as minus the logarithm of Z„, (taken 
along its principal branch defined as log(l) = 0). So 

Y„\ ^ ^ 1 / yp 



A-»--E'o. i-f-hEE, 



p<n '•' V^/ p<„ fc W^ 

These sums clearly converge. Because of this explicit form we can calculate the charac- 
teristic functions. The calculations are done in [181 §3, Ex. 2] and this yields the following. 



ipnit) = E[exp{tt ■ Xn)] = n 2^1 ( ^(«^i + ^2), \itti - t2); 1; - j , 

where t = (^1,^2) G M^, t ■ x = tiXi + ^22^2 is the inner product in M^ and 2^^! denotes the 
Gauss hypergeometric function. Straightforward estimates (see [18] for details) then give 

(1) |v5n(i)| < c(t) exp(— ^(loglogn)|tp), where c is a non- decreasing function (in fact 
one can take a constant); 

(2)¥'n(/^^)^exp(-i|tp). 
The conditions of Theorem H] are fulfilled and hence we have 

for any bounded Jordan measurable set B G C. 

3.9. The Riemann zeta function on the critical line. The results in this section are 
conjectural, but they are of interest to number theorists. By work of Selberg, the central 
limit theorem for logC(l/2 + it) is known, after renormalizing by -^/log log T, see e.g. [9]. 
This is proved by asymptotic estimations of the moments, and there is no known bound for 
the corresponding characteristic functions. Thus, we cannot currently apply our theorems. 
However, Keating and Snaith ([E], [I5]) have proposed the following precise conjecture 
(based on links with Random Matrix Theory) concerning the characteristic function: for any 
t = (ti,t2) G M^, they expect that 

— / exp{it -log ({1/2 + iu))du ^ ^{u)exp( \u\'^] 

T Jo V 4 / 

as T — 7- +00, where the limiting function is the product of the corresponding factors for 
unitary matrices and for the "stochastic" version of (, described in the previous sections, 
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, G(l + ^V^)G'(1 + ^^ 



G(l+2ti) 

p 
(the normalization of logC(l/2 + iw) is obtained by continuation of the argument from the 
value for C(c" + i'^) when a real tends to infinity, except for the countable set of u which 
are ordinates of zeros of C.) 

In [iHl Cor. 9], it is shown that a suitable uniform version of this conjecture implies local 
limit theorems for 

-m{u e [0, T] I log C(l/2 + iu) e B) 

and, as a corollary, implies that the set of values of C(l/2 + iu), m G M, is dense in C, which 
is an old and intriguing conjecture of Ramachandra. 

The vaod-ip framework allows us to show that a much weaker statement than the one 
considered in [18] is already sufficient to get the same local limit theorems. Indeed, we 
consider the following much statement, which of course implies Ramachandra's conjecture, 
as being very likely to be true: 

Conjecture 1 (Quantitative density of values of C(l/2 + it))- For any bounded Jordan- 
measurable subset B d C, we have 



^lim ^ m{u G [0,r] | logC(l/2 + m) G B) 



1 F 
(15) - / exp(it-logC(l/2 + zM))t/n 

^ Jo 



The point is that this follows using Theorem H] from fairly weak decay estimates for the 
characteristic function of logC(l/2 + it) (in comparison with what the Keating-Snaith con- 
jecture suggests). For instance, if for all fc > there exists C^ > such that 

- l + |t|4(loglogT)2 

for all T > 1 and t with |t| < k, then Conjecture [T] is true. Indeed, in TheoremlU we can take 
if to be the characteristic function of a standard complex gaussian and X„ to be (for some 
arbitrary sequence T„ going to +oo) a random variable with law given by the probability 
distribution of logC(l/2 + iu) for u uniform on [0,T„]. These satisfy HI trivially, and H2 
holds with 

An{tut2) = Al{tut2) = yil0gl0gT,(ti,t2), 

because of Selberg's Central Limit Theorem. The hypothesis flTSj) states that, for any A; > 0, 
we have 

\Mt)\<Ckh{Alt) 

for |t| < k, with 

^^'''''^ = TTm"^ 

or equivalently 

MKt)\<Ckh{t) 
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for |S*t| < k. Since h G L^(R^), this gives (|3]), and we get the conjectured statement from 
the local limit theorem. 

The significance of this remark is the fact that, for fixed t 7^ 0, the decay rate of the 
characteristic function which is required is "only" of order (loglogT)"^, which is much 
weaker than what is suggested by the Keating-Snaith conjecture, and therefore might be 
more accessible. 

Appendix A 

We sketch here a proof of Theorem [31 Suppose the support of / is contained in [—k + 
l,k — 1^. Since we can construct approximations for /"*" and /" separately, we can assume 
without loss of generality that / > 0. 

Let e > and let A^ = e~^^^'^'^\ assuming e small enough that N > k. Let then 6' be a 
continuous function on M'^ such that 



e <e <6 + e 



l/(4d) 



and 9 = e + e^/('^'^) on the support of /, while 9 = e outside [—k, k^. Further, let p be 
a trigonometric polynomial in d variables, with periods (2A^, ■ ■ ■ , 2N), which approximates 
f + 9 uniformly on [— A^, N^ up to an error e. Clearly, we have then 

P>f + 9-e>f>0. 

The function p is considered as a periodic function on W^, and it remains non-negative of 
course. The Fourier transform of p, in the sense of distributions, is a finite linear combination 
of Dirac measures, hence has compact support. 

Now we consider 



hix) = n 



sin^faxo 



1 - ^ 
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d 



where a = ^^, for some S > 0; we find that 



ii-sy<h<i 

on [—k, k^. We select 

9{0)-e 



S 



d{9{0)-e+\\f\\^y 

and then we claim that the function gi = ph satisfies gi > /, while J{gi — f)dm can be 
made arbitrarily small if e is small enough. 

Indeed, first of all we have gi > everywhere, while on the support of / we get 

ph-f>{l-6)''{f + 9-e)-f 

>{l-6n9-e)-{l-{l-6r)f 

>il-d6)i9-6)-d6\\f\\^ 

>{9-e)- d6{9 - £ + II / II 00) > 0, by the choice of 6. 

Next the integral J{gi—f)dm is estimated as follows (using the notation C for non-negative 
constants, the value of which may change from line to line, and |a;|oo = maxi<j<d \xj\). On 
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\x\nr, < k, we have 

{gi - f)dm < j (f + e + e- f)dm < {2k)\26 + e^^^'^'^^). 

On the set {k < |x|oo < N}, using the estimate on p, we get 

/ {gi - f)dm < I pdm< 2e{2NY < Ce^''^. 

Finally, where |x|oo > N , we estimate using the bound on h: 
I {gi - f)dm < [ /i(||/|U + e + e)dm 

, + e{x)+e 



l\xi\>N A^a'^+i 



As e — )■ 0, the denominator a'^'^^N tends to oo since 



a 



d+^N = £-1/2 



e 



l/(4rf) \ ('^+l)/2 



and therefore J {gi — f)dm can be made arbitrarily small by choosing e small enough, as 
claimed. 

To conclude, we note that the Fourier transform oi ph is, up to a constant, a convolution 
of the Fourier transforms of p and h. Since the Fourier transform of h is supported on 
[—2a,2aY, the support of the Fourier transform oi ph is therefore contained in the sum of 
two compact sets, which is compact. 

Similarly, using the function f — 6, we construct the approximation function g2 < f. 
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